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Bk from Two-flavor Dynamical Domain Wall Fermions 

Taku Izubuchi ab * for [ RBC collaboration ]t 
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We report preliminary results from an ongoing calculation of Bk for Nj = 2 dynamical QCD with domain wall 
fermions. Simulations have been done with three dynamical quark masses on 16 3 x 32 volumes with L s = 12, 
where the lattice spacing is a -1 = 1.81(6) GeV. Using measurements on ~ 70 lattices for each dynamical mass 
and extrapolating md yn = rn vaj \ to the kaon point, we find B^ s (p = 2GeV) = 0.503(20). 



1. INTRODUCTION 

The kaon B parameter, Br-, is an essential 
bridge between Kq — Kq mixing experiments and 
the CKM matrix. Using domain wall fermions 
(DWF), quenched calculations of Bk have been 
done |1I2I3| . where operator mixing and O(a) er- 
rors are small because of the good chiral proper- 
ties of DWF. Including dynamical fermions is an 
obvious way to improve these calculations. 

Using an improved HMC algorithm for DWF 
1 1I5| . we have generated three ensembles of 
lattices of size 16 3 x 32 for N f = 2 QCD 
with degenerate dynamical masses of mdyn = 
0.02, 0.03, 0.04. Our fermion action is the stan- 
dard DWF action with L s = 12 and M 5 = 1.8 
and our gauge action is DBW2 with (3 = 0.80. 
Extrapolating the rho mass in md yn to the chiral 
limit, we find a lattice spacing of a -1 = 1.81(6) 
and a residual mass m res = 0.00136(5). We mea- 
sure on lattices separated by 50 trajectories and 
have 78, 84, 57 lattices for m dyn = 0.02, 0.03 and 
0.04, respectively 3 . For more details see jS] and 
for results with AS = 1 matrix elements see [H]- 
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2. CALCULATION OF B K 

Bk is defined by 



Br. 



(K°\O LL \K ) 
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(1) 



where A p = sj^d and Oll is the AS = 2 four- 
quark operator (sd)v-A (sd)v-A- For a general 
value of the mass of the pseudoscalars entering 
0, we will use the symbol Bps, reserving Bk for 
the case where the pseudoscalars have the kaon 
mass. 

We first evaluate this by the conventional 
method, where the ratio of a three-point Green's 
function and two pseudoscalar-axial-vector cor- 
relators is taken: 



R(t) 



{J%{h)Q LL (i)j%{h)) 

|(Jf(t )At(t))<At(t)J 5 a (ti)> 



(2) 



(no sum on a). We use Coulomb gauge fixed 
wall sources at to = 4 and t\ = 28 for valence 
masses, m va i, equal to 0.01, 0.02, 0.03, 0.04 and 
0.05. Propagators are found with both periodic 
and anti-periodic temporal boundary conditions 
and averaged to double the temporal length. 

In Figure H w e show R(t). We choose to fit 
data from t = 14 to 17, same fitting range as 
P, to determine B 1 ^. The error bars in this fig- 
ure, and throughout this report, assume that our 
lattices are decorrelated. The fits shown in Fig- 
ure □ give values of Bjfg = 0.557(10), 0.611(6) 
and 0.645(7) for md yn = rn va i — 0.02, 0.03 and 
0.04 respectively. Variations in the plateaus are 
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Figure 1. R(t) for m va i = md yn = 0.02, 0.03 and 
0.04. The lines are the average from t = 14 to 17. 



Figure 2. £?pg for different valence and dynamical 
masses. The filled symbols have m va \ = m.d yn - 



noticeable at the 1-2% level and the assumption 
of 50 trajectories being sufficient for decorrelation 
needs further studies and statistics. However, the 
current data certainly makes a few percent reso- 
lution in the lattice values achievable by possibly 
doubling the statistics. 

In Figure |21 we show our results for £?pg for 
degenerate valence quark masses for each differ- 
ent value of mdyn- The filled points in the figure 
have mdyn = m va i and correspond to the fitted 
lines in Figure ^ (We also have data for non- 
degenerate valence quarks and with our statistics, 
only the average quark mass plays a role.) Our 
simplest result to report for Bk comes from ex- 
trapolating the m va i = mdyn points in Figure [3 
to irif — mj/2 = 0.018. (A pseudoscalar with 
two quarks with m/ = mj/2 has the mass of the 
kaon.) We perform this extrapolation using the 
partially quenched formula 7 

SpsM = ^(l-^(M| s log *jk)) 

+ M^s , (3) 

fitting the three points with m va i = mdyn to 
determine b and bi. We find b = 0.267(12), 
bi = 1.256(79) and = 0.541(9) (statistical 
error only), with a small x 2 /d.o.f. Note that we 
are using the known coefficient of the chiral log- 
arithm, Mpg = 4.06 X (m/ + TOrea),/ = 0.0786 
and A x pt = 1 GeV. The statistical error in the 



lattice spacing gives an additional 2% systematic 
error in the value of B^ , due to the uncertainty 
in the strange quark mass. 

While our data makes extrapolation to m/ = 
m s /2 quite robust, the value of bo is much less 
solid. A small value in the chiral limit has also 
been seen in the quenched calculations, but more 
control over the extrapolation is needed. 

We have also fit the dependence of Bps on 
m va i to @ for each fixed md yn , using the five 
values for m va i. The interpolated values for 
5x(m val = m7/2,m dy „) are 0.537(11), 0.557(9) 
and 0.568(10) for m dy n = 0.02, 0.03 and 0.04 re- 
spectively. We reduce the fitting range to exclude 
the heaviest two points and see the results stay 
same within error. 

An alternative method of extracting Bk is to 
divide the three-point function in the numera- 
tor of QJ by a wall-wall pseudoscalar correlator 
PP. This removes some zero mode effects in the 
quenched approximation, while giving larger sta- 
tistical errors. However here we see no difference 
in the central values for results calculated using 
these two methods. 

3. RENORMALIZATION AND MIXING 

To get to B^ s we need to renormalize B^ 
and control mixing with other dimension six, chi- 
rally disallowed operators. Using NPR JS] in the 
quenched theory we have seen that mixing co- 
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efficients with chirally disallowed operators are 
small. Preliminary NPR calculations on our ~ 40 
dynamical DWF lattices also show that the mix- 
ing coefficients are small (less than a few percent). 

These NPR results are preliminary, since we 
have not fully investigated all the mass depen- 
dence of the procedure and our statistics are 
somewhat limited. However, we can also use the 
f2 parameter introduced in [Q to quantify the chi- 
ral symmetry breaking of DWF and conclude that 
the mixing with chirally disallowed operators is of 
C( TO ros) an d hence negligible. 

We calculate Bk , the renormalization group in- 
variant (RGI) parameter, from B^f via the mul- 
tiplicative renormalization 

B K = Z RGI ' m - MOM (fx) 

x {Z^=\^a)/Z\^a))B^{a) (4) 

A similar equation is used for B^ s except that 

Z RGI,HI-MOM _^ z MS,M-MOM Thc factor 
Z HGI,M-MOM Z ^S=2 ( )/z 2 ( ) ig ca l cu l at ed fol- 

lowing the techniques in jSJ with all quark masses 
equal to md yn - This quantity is extrapolated to 
^dyn = — "ires then the extrapolation to (pa) 2 — 
is done to remove 0(a 2 ) errors. The two loop 
a s (p) formula with N F = 2,A QCD = 300 MeV 
and a -1 = 1.81GcV is used. The uncertainty due 
to the choice of Aqcd is estimated by changing 
Aqcd by ±50 MeV. We check the results stay 
same within error by changing the order of the 
two limits, (pa) 2 — ► and md yn — ► —Tires- 

Following the procedure above, we find 

Z RGI,RI-M0M Z AS=2 /Z 2 j_ 2 g (4 ) and 

z ms,ri-mom z as= 2/Z 2 = 0.93(2) at fi = 2 GeV. 
These numbers agree with the one loop pertur- 
bative calculation, which gives the MS factor to 
be 0.920. 

4. CONCLUSIONS 

Our initial calculation of Bk with dynamical 
domain wall fermions has yielded B^ s (2 GeV) = 
0.503(20) and B K = 0.697(33) . The errors on 
these quantities include the statistical error on 
B^, the error on the Z factors and a 2% error 
reflecting the uncertainty in the kaon mass arising 
from the uncertainty in the lattice scale. These 



Table 1 

Preliminary results for renormalized Bk from va- 
lence quark extrapolations. For m dyn = 0.02, 0.03 
and 0.04 a 2% scale error has been included in the 
error on the renormalized quantities. 



m dyn 


r>lat 
K 


B]f (2GeV) 


B K 


0.02 


0.537(11) 


0.499(22) 


0.692(34) 


0.03 


0.557(9) 


0.518(20) 


0.719(32) 


0.04 


0.568(10) 


0.529(21) 


0.733(34) 


co n 




0.536(6) 




co [21 




0.564(14) 





results correspond to an extrapolation of md yn = 
m va i to m7/2. 

We summarize results for each dynamical lat- 
tice in Table ^ These are obtained by interpo- 
lating m va i to mj/2 at fixed md yn and are some- 
what smaller than quenched DWF calculations at 
a similar lattice spacing. We see a trend in Ta- 
ble U toward lower values as the dynamical quark 
mass is reduced. We have likely underestimated 
our statistical error, due to the long autocorre- 
lation times in dynamical QCD simulations, and 
we have not yet studied the other methods of chi- 
ral extrapolations, the continuum limit or finite 
volume effects. However, it is clear that these dy- 
namical DWF calculations are possible with cur- 
rent computers and hold promise for the next gen- 
eration of computers. 
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